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Abstract. We report results from a search for neutral Higgs bosons decaying to tau pairs 
produced in association with a b-quark in 1.2 fb _1 of data taken from June 2006 to August 
2007 with the DO detector at Fermi National Accelerator Laboratory. The final state includes 
a muon, hadronically decaying tau and jet identified as coming from a 6-quark. We set cross 
section times branching ratio limits on production of such neutral Higgs bosons cj> in the mass 
range from 90 GeV/c 2 to 160 GeV/c 2 . Exclusion limits are set at the 95% Confidence Level for 
several supersymmetric scenarios. 



1. Introduction 

In models with supersymmetry (SUSY), the Higgs sector is expanded relative to the standard 
model (SM). In the minimal supersymmetric standard model (MSSM), two Higgs doublet fields 
are needed instead of one, and the resulting Higgs particle spectrum consists of two neutral 
scalars, a single neutral pseudoscalar and two charged scalars. The ratio of the vacuum 
expectation values of the two doublets is denoted tan/3, and the coupling to weak isospin —1/2 
members of the fermion doublets becomes proportional to tan (3. In the mass range we consider, 
90-160 GeV/c 2 , two of the three neutral higgs bosons will be nearly degenerate in mass and this 
overlap results in a factor of 2 gain in effective cross section. In the MSSM the dominant Higgs 
decay for a pseudoscalar Higgs mass (mi) below «500 GeV/c 2 is to a pair of b quarks, with r 
pair decays occurring with a branching ratio of roughly 10%. However, the (fib — > bbb final state 
suffers from a large multijet background while the (fib — >brr channel offers a much cleaner final 
state, giving the two channels similar sensitivities. Additionally, while the pp — ></> — > rr process 
has a higher cross section than the pp — ► (fib -^brr process, <f> — > rr suffers from irreducible 
Z — ► rr background, while the additional 6-jet in brr largely alleviates this problem, giving brr 
additional sensitivity if ~ Mz ■ 

The DO detector [U EJ [3] has a central-tracking system, hermetic liquid-argon calorimeter, 
and outer muon system. We perform the search in 1.2 fb _1 of data taken from June 2006 to 
August 2007. We search for (fib in the 6rr final state, with one tau decaying to a muon and 
the other to hadrons. This offers the best balance between a large branching ratio and robust 
object identification in the detector. 



2. Event Selection and Background Estimation 

The final state includes one muon, one hadronic tau candidate and one jet coming from 6-quark 
fragmentation, require one isolated muon in the event with px > 12 GeV/c, \rj\ < 2.0 and 
a central track match. Hadronic taus appear as jet-like objects in the detector. We apply 



a neural network (_/ViV T )[3] to distinguish hadronic tau decays from jet fakes. DO has three 
different types of hadronic taus, differing in their number of charged tracks and/or the presence 
of electromagnetic energy in the calorimeter. We make the following kinematic and neural 
network requirements on the hadronic tau candidate: 

• Type 1: E T > 10 GeV/c, p l f k > 7 GeV/c, NN T > 0.9 

• Type 2: E T > 10 GeV/c, ptf k > 5 GeV/c, NN T > 0.9 

• Type 3: E T > 15 GeV/c, 1 track with p T > 5 GeV/c, ^pp > 10 GeV/c, NN T > 0.95. 

Additionally we require at least one jet with pt> 15 GeV/c, \rj\ < 2.5. At least one of the 
jets must be tagged as coming from 6-quark fragmentation by the DO neural network 6-tagger. 
Figure Q] shows the data/background comparison in the visible mass (defined as the invariant 
mass of the muon, hadronic tau, and missing transverse momentum) before and after 6-tagging. 

Backgrounds in this search include VF/Z+jets production, multijet production, Di-boson, 
ttand single top production. The VF/Z+jets contribution is estimated using ALPGEN[5j Monte 
Carlo (MC) interfaced with PYTHIA[6j for hadronization and showering. Signal, Di-boson 
production and tt production are estimated using PYTHIA, and single top production via 
COMPHEP[7l El E]. Multijet background is estimated primarily from data. We use two 
independent methods to estimate multijet production; the first method relies on measuring the 
probability for a jet to be 6-tagged in a multijet-enriched sample, while the second method uses 
the probability for a jet to pass the NN T cut and for a muon to pass the isolation requirement 
in a multijet-enriched sample. We take the average of the two methods as the final multijet 
contribution and include the difference between the two methods as a systematic error. 




Figure 1. Data/background comparison in the visible mass variable (M(/u, r, MET) after muon, 
hadronic tau, and jet selection. Left: before 6-tagging. Right: after 6-tagging. In the legend, 
"QCD" refers to the multijet background. 



3. Multivariate Methods 

After 6-tagging our data sample is dominated by tt and multijet events. We employ two 
multivariate techniques to reject each of the two leading backgrounds. To reject tt background 
we apply a Kinematic Neural Network (KNN) originally developed in [10]. It uses the number of 
jets in the events, the sum of the transverse momenta of the jets (HT), the energy from the four- 
momentum sum of the muon, tau and jets, and the A(ft between the muon and tau candidate as 
input variables. A KNN cut of 0.3 typically offers pa 75% rejection in tt with only a pa 4% signal 
loss. To reject QCD we apply a simple unbinned log-likelihood ratio, trained separately for each 



signal mass point. We consider muon px, tau pr, AR(/j,, r), fj, — r invariant mass, and (/i, r, ^r) 
invariant mass, or visible mass, as input variables, and compute the likelihood of an event to 
be QCD-like or signal-like in all five variables Figure [2] shows the KNN and QCD likelihood 
distributions for signal and background after 6-tagging. 
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Figure 2. KNN output vs. QCD likelihood output with 120 GeV/c 2 signal after 6-tagging for: 
Type 1 taus (left), Type 2 (centre), Type 3 (right.) 



4. Limits and Conclusion 

We cut on both the KNN and QCD likelihood in the 2D distributions of Fig. [21 using expected 
significance as the optimising variable, to determine the cross section limit. Once we have the 
cross section limits, shown in figure O we use FeynHiggs v2.6.2 112 } [T3l 1 14] to interpret the 
limits in four Minimal Supersymmetric Standard Model (MSSM) scenarios: with no m/j-mixing 
and maximal m^-mixing, and with \i = ±200 GeV/c 2 . Fig. U] shows the expected and observed 
limits in the tan/3 vs. uia plane. For uia between 90 and 160 GeV/c 2 we can exclude tan/3 
above 100 in all four scenarios at all mass points. 

We have performed a search for associated neutral Higgs boson production in the brr final 
state with 1.2fb _1 of data with the DO Detector at Fermilab. In the absence of significant signal 
we set cross section times branching ratio limits on <p + b production and interpret the result in 
several benchmark MSSM scenarios. The search is complimentary to bbb searches and provides 
enhanced sensitivity in models where ~ Mz- 




Figure 3. Expected and observed cross section x BR limits vs. Higgs mass. 
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Figure 4. Limits on tan (5 vs. Ma for the no-mixing, /i > case (upper left); no-mixing, /i < 
case (upper right); maximal- mixing, /x > case (lower left); maximal- mixing, fi < case (lower 
right). 
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